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We study experimentally and theoretically polarization-dependent luminescence from an ensemble 
of quantum-dot-like nanostructures with a very large in-plane shape anisotropy (quantum dashes). 
We show that the measured degree of linear polarization of the emitted light increases with the 
excitation power and changes with temperature in a non-trivial way, depending on the excitation 
conditions. Using an approximate model based on the k ■ p theory, we are able to relate this degree 
of polarization to the amount of light hole admixture in the exciton states which, in turn, depends 
on the symmetry of the envelope wave function. Agreement between the measured properties and 
theory is reached under assumption that the ground exciton state in a quantum dash is trapped in 
a confinement fluctuation within the structure and thus localized in a much smaller volume of much 
lower asymmetry than the entire nanostructure. 

PACS numbers: 78.67.Hc, 73.21. La 



I. INTRODUCTION 

Quantum dashes (QDashes) are epitaxially grown 
nanostructures strongly elongated in one of the in-plane 
directions. They can be spontaneously and preferen- 
tially formed in a molecular beam epitaxy process of self- 
assembled growth of, e.g., InAs on InP substrate (Ref. 1 
and references therein). The interest in these structures 
is partly motivated by their applications in telecommu- 
nication lasers due to their higher surface density, fa- 
vorable emission wavelength, wide spectral tunability, 
broad gain, and high speed modulation^—. The exist- 
ing theoretical predictions on some quantum dash prop- 
erties include calculation of their electronic states as well 
as the transition dipole moments and the resulting op- 
tical spectra 5 - - — . Experimentally, the optical properties 
of InAs/InP QDashes have been investigated from the 
point of view of both the ensemble 8 and the single object 
propertie s 9 ' 10 . So far, however, the details of the spa- 
tial character of the confining potential and polarization 
properties of these structures have not been analyzed. 

A scanning electron microscopy (SEM) image of a 
QDash sample, shown in Fig. [TJ reveals shape irregular- 
ities of these structures in the form of sections with lo- 
cally increased thickness as well as zig-zag bends. While 
the shape irregularity is a generic feature of QDash struc- 
tures, the exact morphology varies from sample to sample 
(depending on the details of the growth procedure) and 
either thickness fluctuations or bends can be the dom- 
inating irregularity featureiiii. Both the thick sections 
and the zig-zag corners^, as well as the possible com- 
position inhomogeneities, can act as additional trapping 
centers within the confinement volume of the elongated 
nanostructure. Such an effect of confining in an effec- 
tive potential on length scales smaller than the size of 
the entire nanostructure was observed in V-shaped quan- 
tum wires^. In quantum dashes, the localized charac- 




FIG. 1. Scanning electron microscopy image of a QDash sam- 
ple before overgrowth. 



tcr of emission at low temperatures in low excitation 
regime is suggested by experimental observations on sin- 
gle QDashes, where some features typical for the strong 
confinement regime have been unexpectedly revealed 10 in 
spite of the large volume of these structures. First of all, 
the obtained biexciton binding energy is approximately 
0.5 meV, which is a value characteristic for very small or 
very large quantum objects^. The latter possibility can, 
however, be excluded by noting that the estimated exci- 
ton to biexciton lifetime ratio is close to 2, which is a clear 
fingerprint of a rather strong confinement regim e 15 ' 16 . 

The presence of an additional confinement within a 
larger structure can be seen in the temperature depen- 
dence of the emission lifetime^ or in high resolution 
photoluminescence mapping by using near field optical 



spectroscopy*^. However, these methods require time- 



resolved spectroscopy techniques that suffer from the lack 
of high sensitivity detectors in the 1.55 /im range (nec- 
essary for the investigation of dashes on the InP sub- 
strate) or an AFM-based near field spectroscopy appa- 
ratus. Therefore, a different approach for identifying 
such a two-stage character of the localizing potential in 
QDash structures is highly desired. Such a method can 
be based on the fact that the anisotropy of the confine- 
ment is reflected in the polarization properties of the 
luminescence^—. This relation can be traced back to 
subband mixing effects: Oppositely circularly polarized 
contributions to the optical emission originating from 
heavy and light hole transitions interfere, leading to el- 
liptically polarized emission and to the appearance of the 
preferential orientation of the linear polarizatio n 23 ! 24 . A 
detailed analysis of the carrier states in an elongated 
structure, based on the multi-band k ■ p theor y 25 ' 26 or 
empirical tight-binding approach 2 ^ indeed reproduces the 
polarization properties of the observed emission from 
anisotropic structures. The presumed localizing minima 
in a QDash can be expected to be significantly more 
isotropic than the whole QDash itself. Our idea is, there- 
fore, to use the degree of linear polarization of the con- 
tinuous wave (cw) luminescence of the system as the in- 
dicator of the localization of carriers. 

In this work, we show that indeed the emission from 
the states confined in the trapping centers has much more 
isotropic polarization properties than the radiation origi- 
nating from states delocalized over the whole QDash vol- 
ume due to a considerably different degree of anisotropy 
of their wave functions. For experimental confirmation, 
we investigate the linear-polarization-resolved photolu- 
minescence (PL) from InAs/InP QDashes for different 
carrier distributions obtained by changing the tempera- 
ture and excitation power showing that the polarization 
study of the QDash luminescence can yield information 
on the spatial character of the confining potential and re- 
veal the presence of an additional carrier (exciton) trap- 
ping within a QDash. 

The paper is organized as follows. In Sec. HH we 
describe the system and the experimental setup and 
present preliminary PL results motivating the following 
study. Next, in Sec. MI A[ we present a general the- 
ory relating the polarization properties to the hole sub- 
band mixing, followed by a description of our model of 
a QDash (Sec. IIII B[) . Sec. [TV] presents the results of our 
polarization-dependent PL measurements and theoreti- 
cal modeling. The final Sec. W\ concludes the paper. 



II. SAMPLE AND EXPERIMENTAL SETUP 

The experiment was performed on an ensemble of self- 
assembled InAs quantum dashes grown epitaxially on a 
(001) InP substrate with 3.4% lattice mismatch in the 
molecular beam epitaxy technique. The QDash layer 
is surrounded by additional Ino.53Gao.23Alo.24As qua- 
ternary barriers, lattice matched to InP. The resulting 



nanostructures are significantly elongated in one of the 
in-plane directions, preferentially in [110}^. In Fig. [TJ 
we present an example of a scanning electron microscopy 
image of a QDash sample. It should be kept in mind, 
however, that exact morphology of these structures may 
vary depending on the growth conditions (see Refs. [0 and 
111 ). The lateral dimensions, estimated for a sample as in 
Fig. [JJ are approximately several to tens of nanometers in 
width and even hundreds of nanometers in length, con- 
firming the expected significant shape asymmetry (lateral 
aspect ratio above 4). Their height is typically on the 
order of a few nanometers. This geometry, formed dur- 
ing the self-assembled growth, is a result of a diffusion 
coefficient anisotropy due to the surface reconstruction. 
Further growth details can be found elsewhere 1 - 1 - 3 . For 
the experimental studies, dashes obtained by deposition 
of approximately 1 nm of InAs have been chosen, yielding 
typically the room temperature emission in the range of 
1.55 /im. The surface density of these structures in the 
sample is rather high and exceeds 10 11 cm~ 2 . 

The linear-polarization-resolved photoluminescence 
measurements were performed in a standard photolumi- 
nescence setup. The structures were non-resonantly ex- 
cited by a continuous wave semiconductor laser at the 
wavelength of 660 nm. The excitation powers ranging 
from Pa = 100 nW to 10 4 Po = 1 niW (determined out- 
side the cryostat) were used. The beam was focused onto 
a spot of approximately 0.01 mm 2 on the sample surface. 
The emitted light was dispersed in a 0.5 m focal length 
monochromator and the signal was detected by a liquid 
nitrogen cooled InGaAs CCD linear detector. The PL 
spectra were measured for the two orthogonal polariza- 
tion directions: along and perpendicular to the quantum 
dash elongation axis. From these measurements, the de- 
gree of linear polarization (DOP), defined as 
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h + It ' 



(1) 



was determined, where I\ and I t are the PL intensities 
for the linear polarization parallel and perpendicular to 
the dash elongation axis, respectively. In this work, we 
choose to study the DOP calculated using the intensities 
at the maximum of the ensemble luminescence peak. Us- 
ing intensities integrated over the whole spectral range 
leads to nearly identical results. 

A set of photoluminescence spectra, measured at two 
temperatures for various excitation powers are presented 
in Fig. [21 For each temperature and excitation power, 
two emission spectra for the two orthogonal polarization 
directions ([110] and [110], i.e., perpendicular and paral- 
lel to the quantum dash longer in-plane dimension, re- 
spectively) were measured. The first set of the spectra 
(upper panels in Fig. [5]) have been obtained at low tem- 
perature for three different excitation powers and nor- 
malized, while the second set (lower panels) corresponds 
to a higher temperature. It can be seen that in all the 
cases the polarization along quantum dashes dominates. 
However, at low temperatures, the ratio between the in- 
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FIG. 2. (Color online) Polarization-resolved photolumines- 
cence spectra measured on the ensemble of InAs/InP quan- 
tum dashes for different excitation powers (determined out- 
side the cryostat) and temperatures. The obtained values of 
the DOP at the peak of the PL spectrum are given in the 
figures. The spectra are normalized in such a way that the 
maximum of the signal for the perpendicular polarization is 
the same for all three powers at a given temperature. 



tensities of the linear polarization components increases 
with the excitation power. As a figure of merit we use the 
DOP defined in Eq. ([I}. As can be seen in Fig. ® the ob- 
tained DOP differs by more than a factor of two between 
the lowest and the highest excitation power used. 

Before we proceed to study the polarization-related 
features in detail let us note that, based on the PL spec- 
tra shown in Fig. ® each QDash appears to behave as 
an isolated system, that is, coupling between the dashes 
is negligible. Indeed, if the inter-QDash coupling were 
strong enough to induce considerable carrier transfer on 
the time scales of the exciton life time the emission at low 
excitation powers would be dominated by those dashes in 
which the exciton energy is the lowest, while higher en- 
ergy QDashes would contribute at higher excitation due 
to Pauli blocking (in analogy to the appearance of emis- 
sion from higher shells in a single nanostructure) . This 
would lead to considerable power-dependent broadening 
of the ensemble photoluminescence which is not present 
in the PL spectra (less then 10% increase of the width is 
observed between the lowest and highest powers used in 
Fig.©. 

In addition, only a slight blue shift of the ensemble 
PL feature is seen in Fig. [2] for increasing powers (about 
7 meV between the lowest and the highest powers used) . 
This allows us to exclude heating effects even at the high- 
est powers used, as this would be accompanied with a red 



shift following the band gap reduction according to the 
Varshni law. An additional argument against heating ef- 
fects follows from micro-PL measurements on the same 
sample—, where much more focused beams, hence much 
higher power densities are applied. In spite of this, the 
single QDash lines visible in those experiments show no 
red shift that would be a signature of heating. 



III. THEORY 

In this section, we first develop a general descrip- 
tion of luminescence polarization due to a confinement 
anisotropy in a nanostructure (Sec. lIIl"A"]) and later intro- 
duce a specific model of electron, hole and exciton wave 
functions in a QDash (Sec. IIIIB[) . Although the general 
theory of Sec. IIII Al can be reduced to a simple formula 
that may be able to qualitatively capture the essential 
properties of the system one needs a more accurate char- 
acterization of the actual wave functions, as proposed in 
Sec. IIIIB1 in order to quantitatively account for the ex- 
perimentally observed polarization. 



A. Polarization orientation and hole subband 
mixing 

For the calculation of Coulomb-correlated optically ex- 
cited states of a semiconductor system it is most common 
and convenient to use the hole picture, that is, to de- 
scribe the many-particle states of the almost fully occu- 
pied valence band in terms of the few unoccupied states. 
Single particle (electron and hole) states confined in a 
nanostructure can then be approximately described in 
the single-band envelope function approximation by the 
envelope wave functions tp c ( v )i\(r) , where A denotes the 
subband, i labels different eigenfunctions of the confine- 
ment potential, and c(v) refers to the conduction (va- 
lence) band. Note that in this notation ip V i\ is the hole 
wave function (rather than an electron wave function for 
a valence bad state). On the other hand, for multiple- 
band calculations, like the k ■ p-based perturbation the- 
ory to be presented below, as well as for a formal dis- 
cussion of inter-band transitions, it is by far more con- 
venient to use the electron picture, in which the states 
of electrons in the valence band are directly represented 
by their wave functions. The hole wave function and the 
valence band electron wave function are simply related 
by complex conjugation. Thus, in the electron picture, 
the envelope wave functions for the valence band elec- 
tron states are ij}* iX {r). The corresponding creation and 
annihilation operators (in the electron picture) will be 

denoted by o; (v)iA , a c ( v ) iX . 

If the Coulomb interaction and hole subband mixing is 
included the energy eigenstates of a confined exciton can 



be written in the general form 
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where e x ,e y are the unit vectors along the crystallo- 
graphic axes. From Eqs. ©, §5§ and ([5]), the corre- 
sponding components of the interband dipole moment 
are 



where A and A' run through the valence and conduc- 
tion subbands, respectively, j3 indicates different exci- 



ton states, c 
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are the coefficients for the expan- 



sion in single-subband single-particle states, and |0) is 
the ground state of the crystal. 

The polarization of the light emitted in the recombi- 
nation process from a conduction band state to a valence 
band state is determined by the interband matrix element 
of the positive frequency part of the dipole moment op- 
erator d (Ref. 1281) . 
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d = d\x> j d 3 rip wi \{r)ip c j\i(r)a\ iX a c j\<, 

where d\\> is the interband matrix element between the 
states at the T point of the Brillouin zone in the bulk 
material (note that, according to our notation, ipvi\{f) 
is the complex-conjugated valence band wave function). 
The dipole moment for the optical transition is therefore 
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Let us assume that the initial state is spin-up (A' = 
+1/2). The only non- vanishing in-plane components of 
the dipole moment are to the +3/2 heavy hole (hh) band 
and to the —1/2 light hole (lh) band, with the corre- 
sponding bulk matrix elements along the crystallographic 
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Note that the magnitude of the dipole matrix elements, 
do, is irrelevant for the present study since the DOP is a 
relative quantity, as follows from Eq. flj . 

Since the studied system consists of QDashes elongated 
in the [110] direction we define unit vectors parallel and 
transverse to this direction, 
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If the average occupation of an exciton state with en- 
ergy E is n(E) then the intensities of light emitted from 
the state /?, polarized parallel and perpendicular to the 
direction of the QDash elongation are proportional to 
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tively 



respec- 



In order to provide a description relevant to an in- 
homogeneously broadened ensemble we consider a set 
of QDashes with variable parameters that we formally 
jointly denote by a single symbol 77, representing the full 
set of relevant parameters. The random distribution of 
the QDash parameters is represented by the distribution 
function f{rj). At this point, the number and the physical 
nature of the variable parameters as well as the distribu- 
tion function can be arbitrary but a simple choice will 
be proposed in Sec. IIIIBI in order to perform numerical 
calculations. 

For the PL response recorded from an inhomogeneous 
QDash ensemble at a single emission energy E, the PL 
intensities at the two linear polarizations are then 



h/t(E) = 
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where E v p is the energy of the state \X^) and d[? t are 
the corresponding components of the dipole moment in 
a QDash characterized by the parameters r\. According 
to Eq. fT) and using Eqs. © and ©, the ensemble DOP 
at the emission energy E is then 
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where ct^y 1S ^ ne oscillator strength parameter, calcu- 
lated according to Eq. ([4]) for a QDash with parameters 
r\ and n^iE^fi) are the average occupations of the en- 
ergy levels in such a QDash (occupation redistribution is 
assumed to be independent for each QDash, that is, no 
redistribution between the dashes takes place). As can 
be seen from the above equations, mixing between light 
and heavy hole states leads to preferential polarization 
of the emitted light, which depends on the relative phase 
of the two contribution a 23 ' 24 . 

In our theoretical modeling, we assume that the lowest 
hole state consists mostly of a heavy hole component with 
some admixture from light holes, as is typical for self as- 
sembled structures, where the strain is compressive. This 
is justified in many systems since the light hole states are 
shifted in energy with respect to the heavy hole states 
due to confinement and strain, while the inter-subband 
coupling elements are relatively small. In order to derive 
the light hole contribution one uses a perturbation theory 
based on the k ■ p (Kane) Hamiltonian 3 ^. As can be seen 
from the structure of the Kane Hamiltonian, the heavy 
hole +3/2 state is coupled in the leading order to both 
light hole states. However, the exciton state involving 
a spin +1/2 electron and a spin +1/2 hole is dark and 
does not contribute to the optical properties in the lead- 
ing order. The calculations of the degree of polarization 
[Eq. (|8j)] have been performed with an assumption that 
coupling to the +1/2 light hole states can be neglected. 

The part of the Kane Hamiltonian relevant to the light 
hole admixture is 



V = 
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where R = i?k + Ra is the matrix element of the Kane 
Hamiltonian coupling the spin 3/2 heavy hole subband 
with the —1/2 light hole subban d 25 ' 30 . The kinetic con- 
tribution to the inter-subband coupling is 
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where kj = —id/dxj, too is free electron mass and jj are 
the Luttinger parameters. The strain-induced contribu- 
tion is 
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where b v and d v are valence band deformation po- 
tentials and £ij are the strain tensor components. 
Note that we use the standard definition of the basis 
functions^ -, consistent with the general theory of the an- 
gular momentum 3 ^,, which differs from that used in many 
papers employing the k ■ p theor y 25 ' 32 . 

Assuming the zeroth-order exciton state to be of purely 
heavy hole type with, say, a spin-up electron, 
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the first order correction due to the inter-subband cou- 
pling is approximated as 
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where A£ih is the average energy separation between 
heavy and light holes. Substituting this to Eq. (|4|) one 
finds 
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are the kinetic and strain-induced contri- 
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for x = k, s, where we used the completeness relation for 
the wave functions ipvk, -i/2( r )- 



As the wave functions can be chosen real, a 
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is real and only the imaginary part of «_, ,„ i/ 2 > re ~ 
lated to the imaginary (anti-hermitian) parts of i?k and 
R s , contributes to Eq. ([8]). In the leading order in the 
subband-mixing term one then finds the DOP in the 
form n = Hk + n s , where we split the effect into the ki- 
netic and strain-induced contributions defined by Eq. ((SJ) 
with OL_i/n I/O replaced in the numerator by its kinetic or 

strain-related parts &_i/ 2 i/ 2 , % = k,s, respectively, and 
neglected in the denominator. 

The imaginary (anti-hermitian) part of i?k can be writ- 
ten in the coordinate frame related to the QDash elon- 
gation as 
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where k\ — (k x — k y )/\/2 and fc t = (k x + k y )/\/2. From 
Eqs. © and © it is then clear that the DOP is de- 
termined by the matrix elements of kf and fc 2 between 
the states involved in the optical transition, hence by the 
symmetry of the wave function. 

In order to estimate the strain-related part one would 
need to know the exact strain distribution in the nanos- 
tructure, which is beyond the scope of the present paper. 



Nonetheless, if the localization is due to relatively small 
and smooth shape fluctuations then the strain field can 
be expected to vary only very weakly over the volume 
of the QDash (note that the situation can be completely 
different if the trapping takes place in corners or bends 
on the QDash where the strain distribution is likely to 
be much different from that in the straight segments). 
Here, we will adopt the assumption of a nearly constant 
strain field. Then, comparing Eq. ([9]) for R s — const with 
Eq. (T5J) one immediately gets 
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The strain-induced contribution to the DOP in this case 
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and is the same for every state. 

While the result contained in Eqs. |(5J| and §§§ is rather 
involved and requires a numerical solution for the exci- 
ton states, the essential features of our theoretical predic- 
tions are already clear in a simplified, semi-quantitative 
treatment^. If the Coulomb coupling is neglected then 
for a single bright exciton state with both the elec- 
tron and the hole in the single particle state «o one has 
c (i 3/2)(i 1/2) = buo&jia- in this case, if the electron and 
hole envelope wave functions are assumed identical then 

4/2,1/2 = 1 and 
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where 



denotes the average value in the envelope 



state io- Typically (assuming approximately box- like 
confinement model), 
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where L and D are the confinement dimensions along 
and perpendicular to the QDash, respectively and n\, n t 
are the quantum numbers associated with the excita- 
tions along and perpendicular to the QDash elongation. 
Therefore, from Eq. ©, one can estimate the kinetic con- 
tribution to the DOP of light emitted from a single state 
in this case 
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Of course, for other confinement models, the quantum 
numbers will enter in Eq. (|ll|l in a different way (e.g., 
linear for a harmonic oscillator). The main conclusions 
will, nonetheless, be the same. 

The approximate formula (|Iip captures the essential 
effect of the confinement shape on the polarization prop- 
erties of the emitted light. For an in- plane isotropic con- 
finement (L = D), there is no preferred polarization axis. 



For an anisotropic structure, the light is preferentially 
polarized along the structure, to the degree dependent 
on the occupation of the higher excited states. It turns 
out that this simple model is sufficient to account for 
the high-temperature spectrally resolved degree of polar- 
ization for an ensemble of QDashes^. However, for the 
study of more subtle effects, like additional localization, 
one needs a more detailed approach, including a real- 
istic confinement shape, the Coulomb correlations, and 
the thermal distribution of occupations. Therefore, in 
the following section we introduce a more detailed and 
accurate model of electron and hole confinement as well 
as Coulomb-coupled exciton states. While Eq. (fTTj) will 
be useful for the qualitative interpretation of the results, 
the general Eq. J8]) with the wave functions derived in 
Sec. MI Bl will be used in the numerical modeling of the 
polarized emission that will be compared with the exper- 
imental results. 



B. Model of a QDash 

In order to obtain quantitative estimates of the de- 
gree of polarization of the radiation emitted by a QDash 
structure occupied by a single exciton we now introduce 
a simple model of the system. This will allow us to cal- 
culate the kinetic contribution to the polarization, which 
critically depends on the character of the wave functions. 
The strain-related part n s is assumed to be similar for 
all the states as discussed above. Therefore, this part 
will be treated as a constant parameter to be found from 
fitting to experimental results. 

Our choice here is to model the carrier trapping 
by width fluctuations, which are the dominating type 
of confinement irregularity in some InAs/InP QDash 
samples 11 . In the Appendix, we show that localization 
appears also on the zigzag corners. In fact, the exact lo- 
calization mechanism is not essential for our conclusions. 
Our choice is convenient for our present purpose as the lo- 
calization on geometrical fluctuations can be consistently 
modeled within a simple approach, while trapping on the 
bends may involve nontrivial strain fields and, therefore, 
require much more detailed structural modeling, which 
is beyond the scope of this paper. 

The Hamiltonian for a single carrier (electron or hole) 
confined in a QDash in the single-band effective mass and 
envelope function approximations is 



h 2 

2m* 



V(r), 



where a denotes the type of a carrier (a = e, h) . The 
QDash confinement potential V(r) is modeled as a 3D 
potential well with the shape reproducing the essential 
features of the QDash geometry, in particular the pres- 
ence of a widening that can trap the carriers (see Fig. [3]). 
We assume that the QDash structures have the cross- 
section in the form of a circular segment, with the base 
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FIG. 3. (a,b) A schematic representation of the QDash geom- 
etry used in our model: top (a) and section (b) view, (c) An 
exact plot of the QDash width as a function of position along 
its length (only the central section of the structure is shown) . 



width (chord length) changing along the QDash length 
(the x coordinate) according to 



D{x) = D a 



AD(l + 4e 



-b\ 



l + ie-bcoshfibx/xo)' 



where Dq is the QDash base width away from the widen- 
ing, AD is the depth of the fluctuation, xq is the length 
of the fluctuation, and b defines the shape of the widening 
(we choose b = 20; see Fig.^c) for the shape of D(x)). fn 
our model, the widening of the structure is located sym- 
metrically in the center of the QDash (preserving the D 2 
symmetry of the structure) . A non-central location of the 
widening modifies the details of the spectrum due to dif- 
ferent selection rules but does not lead to essential modifi- 
cations of the polarization properties discussed here. The 
height of the dash is H(x) — (D(x), where £ is a constant 
height to width ratio. For the numerical computations we 
use standard material parameters for the InAs/InP ma- 
terial pair— and parameter values computed for strained 
nanostructures in this material system 3 ^: The effective 
band offsets between the materials (including strain ef- 
fects), defining the depth of the confinement potentials, 
are taken as 400 meV and 250 meV for electrons and 
holes, respectively. The effective masses used in our cal- 
culations are 0.037mn for electrons and 0.33mn for holes, 
where mo is the free electron mass. The Luttinger pa- 
rameter is 73 = 9.29. 

The approximate single-particle envelope wave func- 
tions for electrons and holes are found using a variational 
method generalizing the "adiabatic" approximation 36 , 
using the fact that the confinement along the x direc- 
tion is much weaker than in the other two directions and 
changes smoothly. First, for each x, we variationally min- 
imize the single-particle Hamiltonian (where we suppress 
the carrier type index a for clarity) 
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d 2 d 2 



2m* \dy 2 dz 2 / 
in the class of 2-dimensional harmonic oscillator (2DHO) 
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FIG. 4. (Color online) The one-dimensional effective potential 
along the QDash eo(x) and the spectrum of single-particle 
levels for the hole (a) and electron (b). 



ground state wave functions 



i>o{y,z;x) = 

1 

y/l z (x)ly(x)w 



exp 



z {x)f 



2l 2 (x) 



2P y (x) 



(12) 



with the variational parameters l y {x) and l z (x) corre- 
sponding to the wave function widths in the y and z 
direction and zq(x) representing the position of the cen- 
ter of the wave function along z. Next, we generate the 
set of x-dependent effective potentials 



dz / dyip*(y,z;x)H yz ip n (y,z;x), (13) 



.(*) 



where ip n (y, z; x) is the wave function of the 2DHO (with 
the same parameters as those obtained variationally for 
the ground state) representing the nth state along y. 
This allows us to approximately account for the excited 
states along y (however, we restrict the dynamics in the 
strongest confinement direction z to the ground state). 
The approximate energies e n {x) are then treated as effec- 
tive potentials for the one-dimensional eigenvalue equa- 
tions in the QDash elongation direction, 



d 2 



2m* dx 2 



e n {x) 



Jnm[X) — *- J nmJnm\%)j V J 



which are solved numerically (separately for each n) . The 
full envelope wave functions are then 



i>nm(x, y, z) = lj) n (y, Z] x)f nm (x) 



(15) 



In the following, both quantum numbers n and m will be 
denoted by a single index i. 

In Fig.|U we show the one-dimensional effective poten- 
tial along the QDash en (a;), corresponding to the lowest 
states in the yz plane, and the spectra of single-particle 
electron and hole levels. As can be seen, for the param- 
eter values chosen here, there is one trapped state both 
for the electron and for the hole. 



Based on the calculated single-carrier states, one can 
construct the product basis for the excitonic states and 
diagonalize the system described by the Hamiltonian 



i i 

+ ^2 V ijkia\h f J h k ai, 

ijkl 



(16) 



where aj, a^ and h\, hi are electron and hole creation and 
annihilation operators, respectively, E^' are the energies 
found from Eq. {T3J) for electrons and holes, and Vijki are 
the matrix elements of the electron-hole interaction, 



V; 



:-h 

ijkl 



(ij\H e _ h \kl) 



- 1 



An esq \r e -r h 



■?Pck(rh.)ipvi(r e 



where Eq is the vacuum permittivity and e is the rela- 
tive dielectric constant of the QDash material (e = 14.6 
for InAs). Upon the diagonalization of the Hamiltonian 

given by Eq. (fT6|) one obtains the coefficients c/* 3 , 2 \ia 1 i 2 \ 
defining the exciton eigenstates which are then used in 
the calculation of the DOP according to Eq. © (in the 
leading order in subband mixing). 

The inhomogeneous ensemble of QDashes is modeled 
by assuming a Gaussian distribution of the sizes (cor- 
responding to the formal distribution function f(rj) in- 
troduced in Sec. IIII A[) with the standard deviation of 
10% and with a proportional scaling of all the dimen- 
sions. The average length of a QDash is taken to be 
L = 150 nm, the QDash lateral aspect ratio is L/D = 6, 
the height to width ratio is £ = 1/5.5, and the size of the 
width fluctuation is xq = 8 nm, which corresponds to a 
typical geometry of real structures*^". In our numerical 
modeling, AD /Do and AE\h are adjustable parameters. 

For the numerical calculations, we choose a set of sin- 
gle carrier states with lowest eigenenergies. In order to 
achieve convergence we restrict our single carrier states 
to about 25 electronic states (with the dynamics along 
the transverse direction restricted to the ground state), 
and 50 hole states. Since the effective mass of heavy 
holes is larger than that of an electron, also states with 
excitations along y direction (n = 2), which are coupled 
to the n = states by the Coulomb interaction, are in- 
cluded. Depending on the exact values of the QDash 
shape parameters, the number of states included in the 
calculations may slightly vary. Luminescence at finite 
temperatures is calculated using Boltzmann distribution 
for the occupation of exciton states. This is equivalent 
to the assumption that thcrmalization is fast enough to 
assure full relaxation to equilibrium at the lattice tem- 
perature on time scales much shorter than the exciton life 
time. This should be the case in the system under discus- 
sion since the spectrum of a strongly elongated structure 
is rather dense and the optical response originates from 
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FIG. 5. (Color online) (a) Dipole moments for optical transi- 
tions for the exciton states obtained from numerical calcula- 
tions for a 150 nm long QDash with D/ Do = 0.1. (b) and (c) 
Emission spectrum from a single QDash taking into account 
thermal distribution of the carriers at two different tempera- 
tures. Here, relative intensities (normalized to the intensity 
of the ground state transition) are shown. 



structures where electron-hole pairs are captured, which 
enables the efficient electron-hole scattering-assisted re- 
laxation channels. 

In Fig. O we show the optical properties of a sin- 
gle QDash obtained from our numerical calculations. 
Fig. [5ja) presents the interband dipole moments (oscilla- 
tor strengths) for various exciton states confined in the 
QDash. In Fig. OJb,c), the computed emission spectrum 
at finite temperatures is shown. The lowest energy line 
corresponds to a state trapped in the width fluctuation 
while the higher energy lines are related to states delocal- 
ized in the whole QDash volume. The lowest delocalized 
states have larger oscillator strengths than the trapped 
ground state but because the latter is down-shifted by 
the trapping and Coulomb energy the contribution from 
the delocalized states becomes large only at relatively 
high temperatures. Note, however, that the appearance 
of the higher energy emission lines in a single QDash 
spectrum some 10 meV above the fundamental transition 
does not considerably affect the ensemble emission shown 
in Fig. [2] because of the large inhomogeneous broadening 
of about 50 meV. As a result, the ensemble emission fea- 
ture is rather symmetric at all temperatures, reflecting 
the symmetric ensemble distribution of QDash morphol- 
ogy features, which is common for this material system^. 



IV. RESULTS AND DISCUSSION 

In this section we present the results of measurements 
and theoretical modeling that aim at a complete charac- 
terization and explanation of the polarization properties 
of the luminescence from QDashes, as a function of tem- 
perature and excitation power. It should be noted that 
our theory presented in Sec. IHII is formally limited to 
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FIG. 6. (Color online) Points: Temperature dependence of 
emission degree of polarization from the ensemble of InAs/InP 
quantum dashes in three different excitation regimes as indi- 
cated. For clarity, the error bars are shown only for the lowest 
power data where the measurement uncertainty is the largest. 
Solid line: the theoretical result. 



the low power regime, where a single exciton model for 
a QDash emission is valid. Taking into account a rea- 
sonable estimate of losses, the power density and laser 
spot size as given in Sec. [TT] corresponds to the photon 
fluence rate on the order of 10 6 /(cm 2 -ns). With the av- 
erage exciton lifetime of 1 ns, one gets an estimate of the 
QDash occupation on the order of 10~ 5 (assuming the 
exciton lifetime of 1 ns) for the lowest excitation power 
used. Double occupation under weak excitation condi- 
tions is therefore highly improbable, even if one allows 
for the extension of charge life time due to dark spin 
configurations and trapping of unpaired carriers. Hence, 
the single-exciton theory provides a basis for the inter- 
pretation of the temperature dependence of the DOP at 
low powers. However, qualitative or semi-quantitative in- 
terpretation of power-dependent spectra is also possible, 
based on the concept of state filling (due to Pauli block- 
ing) and on the general tendencies predicted by Eq. (TJTJ). 

In order to systematically study the dependence of the 
DOP on the excitation conditions and temperature we 
have measured the DOP as a function of temperature 
for three different excitation powers, shown in Fig. [5] 
When a very low excitation power is used the DOP at 
low temperatures is on the level of 0.12-0.13 and then 
increases with temperature up to about 0.27 at above 
60 K (which roughly corresponds to 5 meV activation 
energy). The dependence is similar for the intermediate 
excitation power but the initial (low temperature) value 
is higher (about 0.19). Eventually, the DOP is almost 
constant on the level of 0.27 over the entire temperature 
range for a very high excitation. 

These results are consistent with the emission prop- 
erties from either locally trapped or delocalized (in the 
whole QDash volume) states. Indeed, as discussed in 
Sec. IIII1 the DOP for the emission from a given state 



depends on the imaginary part of the relevant matrix el- 
ement of the fc • p Hamiltonian given by Eq. (p~0|) which, 
in turn, is determined by the values of the wave vec- 
tor components k\,k t (along and perpendicular to the 
QDash elongation) characteristic for the relevant wave 
function. At low temperature and low excitation, the de- 
tected emission originates mainly from the ground state 
excitons localized at the QDash potential fluctuations, 
which are more isotropic than the dash itself. Due to the 
reduced anisotropy, k\ ~ k t , hence the kinetic compo- 
nent to the DOP is small. This results in a small DOP, 
similar to that observed for slightly asymmetric quantum 
dots^I. By increasing the temperature, the excitons are 
thermally released into the whole QDash volume. For 
such low-energy extended states the confinement along 
the QDash is much weaker than across the structure. 
Since k\ ~ 1/L and kt ~ l/-Do these two wave vector 
components become imbalanced and hence the emission 
becomes more linearly polarized. This effect is also clear 
from the simplified Eq. (TTTJ) if one substitutes for L the 
effective confinement size for a given state (L ~ D for 
a trapped ground state, L 3> D for a delocalized higher 
energy state). 

With increasing excitation power, the higher states 
with anisotropic wave functions become filled and con- 
tribute to luminescence already at low temperatures. 
Therefore, the low temperature DOP increases and the 
amplitude of the temperature-dependent change is re- 
duced until, for the highest excitation, the contribution 
of the localized excitons to the emission becomes neg- 
ligible and the DOP does not significantly change with 
temperature. Remarkably, the saturation level for the 
DOP is the same for all the excitation powers, i.e., it 
is insensitive to the excitation conditions, reflecting the 
intrinsic asymmetry of the carrier states trapped in the 
QDash. 

The results obtained from the full theoretical model 
presented in Sec. IIIII using the exciton wave functions 
found in Sec. IIIIBI and the general equations © and 
© for a certain set of parameters are shown by the 
green solid line in Fig. [5] As our theory is restricted 
to single-exciton states the theoretical curve corresponds 
to the low excitation limit. The S-shaped temperature 
dependence is reproduced if one assumes a QDash en- 
semble with an appropriately chosen relative amplitude 
of the widening. For the curve in Fig. [5J AD /Do =0.1. 
For this value of the shape parameter, only the exciton 
ground state exhibits highly localized properties, which 
are reflected by a moderate value of the calculated DOP 
for low temperatures. Clearly, the model is able to re- 
produce all the qualitative features of the measured de- 
pendence: When only the ground state is occupied the 
degree of polarization of emitted radiation is small. At 
higher temperatures, emission from higher energy exci- 
tonic states (extended over the whole area of a QDash) 
contributes, leading to the increased value of DOP. The 
theoretical result shown in Fig. [5J is obtained for the sep- 
aration between light end heavy holes AE^ = 20 meV 
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FIG. 7. (Color online) Influence of the amplitude of the 
widening of a QDash on the calculated DOP of the system. 
The strain contribution is assumed the same in all the cases. 



(the amplitude of the DOP change scales inversely pro- 
portional to this parameter). While this value is rather 
small, one should notice that this energy is an effective 
parameter averaged over all the heavy hole states con- 
tributing to the luminescence and should not be treated 
as the separation between the lowest heavy hole state and 
the light hole subband (in the decoupled band picture). 

As can be seen, our model is also able to account to 
some extent for the observed decrease of the DOP at 
high temperatures for lower excitation powers. This ef- 
fect can again be understood with the help of the simpli- 
fied Eq. I|lip. As the temperature grows, the first process 
is the release of carriers from the additional confinement 
center which leads to an increase of the effective confine- 
ment length L in Eq. ([TTj) as discussed above. At first, 
the lowest lying states are filled so that m ~ 1. However, 
when the temperature further increases higher excited 
states are occupied with growing values of n\, while rit 
remains small due to the stronger confinement and higher 
excitation energy in the transverse lateral direction. It is 
clear from Eq. (TPTJ) that this increase of n\ reduces the 
value of the DOP as indeed observed in the low-power 
experimental data and in the theoretical result. Quanti- 
tatively, the strong decrease of the low-excitation DOP at 
high temperatures is not reproduced by the theory prob- 
ably due to the simplified nature of our QDash model 
which does not properly account for the states in the 
higher energy sector. Moreover, one needs to keep in 
mind that the signal values at these conditions are very 
low, hence the uncertainty of the measured DOP is very 
high in this range (see error bars in Fig. [5]). 

Let us emphasize that the additional localization is 
essential for the explanation of the observed polariza- 
tion properties of the PL from QDash structures. The 
experimental results cannot be reproduced even qualita- 
tively in the absence of additional confinement inside a 
QDash. When the shape of a QDash is more uniform the 
ground state of the system is not trapped and has a rather 



FIG. 8. (Color online) The comparison between the temper- 
ature (bottom horizontal axis, blue triangles) and excitation 
power (top horizontal axis, red circles) dependence of emission 
degree of polarization for the ensemble of InAs/InP quantum 
dashes. 



anisotropic wave function, similar to all the other states 
in the structure. As shown in Fig. this results in a 
relatively high degree of polarization already at low tem- 
peratures, which decreases at higher temperatures due 
to increasing occupation of higher excited states. This 
kind of behavior is a fingerprint of a nanostructure with 
a uniform width. When the widening of a QDash is 
large enough to confine carriers the results change con- 
siderably. Thus, the character of polarization properties 
of emitted radiation reflects the shape of confining po- 
tential: Different confinement conditions lead to quali- 
tatively different polarization properties of luminescence 
from the system. Interestingly, the saturation level and 
the high-temperature behavior are different for different 
sizes of the shape fluctuation, which shows that the pres- 
ence of the trapping potential is important even when the 
emission from the delocalized states becomes dominant. 
Effects similar to those observed in the temperature 
dependence can be seen in the excitation power depen- 
dence of the linear polarization (red circles in Fig. [H]). 
The DOP increases and finally saturates on the same 
level as in the temperature dependence. Obviously, the 
cases of high temperature and high excitation power are 
physically different: In the former, we deal with a single 
exciton system in which the occupation of higher states 
is due to thermal redistribution. In the latter, one has a 
many-body system with excited state occupation forced 
by state filling. The correspondence between the DOP 
observed in these two essentially different cases provides 
a strong confirmation that the effect is due to increased 
occupation of excited states, which is a common factor 
of these two physical situations: Since the higher states 
have anisotropic wave functions they contribute strongly 
to polarized emission. A straightforward comparison of 
temperature and excitation power influence on the DOP 
is presented in Fig. [5] where the DOP axis is common 
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and the temperature and excitation power axes (horizon- 
tal ones) are translated with respect to each other and 
scaled in order to emphasize the correspondence between 
the shapes of the curves. This figure shows directly that 
the two external factors give a similar effect if the exci- 
tation power values are presented in a logarithmic scale. 
This is consistent with the fact that the occupations of 
excited states scale with the excitation power according 
to a power law, while their dependence on temperature 
is exponential. Therefore, the agreement between the 
two dependencies supports the idea that the distribution 
of carriers over the confined states in the QDash is cru- 
cial for the polarization properties of the emitted light: 
Higher excitation power creates more electron-hole pairs 
which occupy higher energy states according to Pauli ex- 
clusion principle, whereas elevated temperatures result in 
the redistribution of carriers already existing in the sys- 
tem allowing them to occupy higher energy levels. Es- 
sentially, however, both mechanisms lead to the observa- 
tion of radiative recombination of carriers from higher- 
energy states that are expected to be confined in the 
whole QDash volume. At the same time, the logarith- 
mic scaling between the temperature and power depen- 
dence is a strong argument in favor of state filling as the 
mechanism of the power dependence, as opposed to heat- 
ing. If the excitation power acted via heating the local 
temperature would be proportional to the power and the 
dependencies on the power and temperature would scale 
linearly, which is not the case. 



operation of some optoelectronic devices, especially in 
those applications where polarization control is impor- 
tant (e.g., polarization- insensitive optical amplifiers^). 
The strongly confined character of the lowest carrier 
states in QDashes will also have a considerable impact 
on the performance of futuristic photonic devices which 
are based on quantum electrodynamics experiments^ - — . 
It will affect, for instance, the exciton fine structure and 
anisotropy splitting energies, the energies and even rel- 
ative spectral position of the excitonic complexes in a 
QDash, and finally the kinetics of the transitions be- 
tween the excitonic states. The effect of carrier (exci- 
ton) trapping will affect both the total exciton oscillator 
strength and the polarization selectivity which become 
crucial when QDashes are used as quasi-zero-dimensional 
emitters placed inside a microcavity, suitable for sin- 
gle photon sources in the telecommunication wavelength 
range of 1.3-1.55 /im. 
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V. CONCLUSIONS 

We have investigated polarization-resolved lumines- 
cence from an ensemble of InAs/InP quantum dashes. 
We have developed a theory relating the observed degree 
of linear polarization to the asymmetry of the carrier 
wave functions via the anisotropy effect on the hole sub- 
band mixing. By comparing the temperature and power 
dependent degree of polarization to the theoretical pre- 
dictions we were able to conclude that the lowest carrier 
states usually have a much lower degree of asymmetry 
than the structure itself, which can be related to addi- 
tional trapping of the excitons to potential fluctuations 
within the dash volume. Thus, our findings reveal a non- 
trivial character of carrier states in these systems of tech- 
nological and applicational relevance. 

Our analysis shows, in addition, that the cw 
polarization-resolved spectroscopy can be used as a probe 
of the localization effects within such elongated objects. 
Thus, this widely available optical tool can yield impor- 
tant information on the actual properties of carrier states 
in anisotropic systems. 

Apart from the general interest in the electronic 
structure of quantum dashes, the knowledge on their 
polarization properties and the nature of the carrier 
states is essential for modeling the emission properties 
of such nanostructures and can be important for the 



Appendix A: Localization on zig-zag corners 

In this Appendix, we verify that zig-zag bends act as 
localization centers in a similar way to local widenings. 
Thus, the theory presented in this paper is applicable also 
to QDash systems where bends, rather than widenings, 
are the dominating kind of shape irregularities. Here, we 
aim at demonstrating the localization effect qualitatively, 
therefore we present a simple model, including only the 
geometrical features of the confinement. It should be 
kept in mind, however, that in a real system a specific 
strain distribution will occur in the bend area, which 
must be taken into account in more realistic modeling. 

We consider a QDash of dimensions as in the main 
body of the paper with a homogeneous cross section but 
with a 90° bend present in the middle (see the outlines 
in the upper part of Fig. [9]). The electron confinement 
is modeled as a three-dimensional potential box with the 
confinement potential depth and cross-section geometry 
as described in Sec. IIIIBI The electron wave functions 
are found in a way similar to that presented in Sec. IIIIBI 
(see also Fig. [3]) but this time only the one-dimensional 
Schrodinger equation along the strongest confinement di- 
rection z is solved, yielding a two-dimensional effective 
potential for the eigenvalue problem in the xy plane. The 
latter is solved by expansion into plane waves. 



12 




50 nm 




In Fig.[9la), we show the resulting probability densities 
for the ground and first excited electron state for the bent 
dash, as compared with a straight, rectangular dash. It is 
clear that the ground state gets localized on the corner, 
while the first excited state remains delocalized in the 
same way for both the bent and straight QDash. 



ground state 



1 st excited state 



FIG. 9. (Color online) The electron densities in the xy plane 
(integrated over z) for the ground and first excited state of a 
bent and straight dash. 
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